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Abstract 

It is shown that the commonly accepted relationship between the Landau singularity 
in the running coupling constant of QED or QCD and the renormalon singularities in the 
Borel sums of perturbation theory expansions is only a particular feature of the restriction 
of the perturbative /3-function to the one loop level. 
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1. The success of the Standard Model in describing particle physics at present energies 
relies heavily on the perturbation theory approach as applied to the underlying gauge field 
theories : SU(S) x SU(2)l x U{1). Yet, it is known since the early work of Dyson Q that 
perturbation theory itself suffers from ambiguities which originate in the fact that physical 
quantities are not analytic in the coupling constants which define the expansion parameters. 
This is reflected in a characteristic /c!-growth pattern in the coefficients having the same 
sign of the large order terms in the perturbation series []. In theories like QED this growth 
originates in the large momentum integration region of virtual photons dressed with vacuum 
polarization corrections and leads to singularities in the associated Borel plane; the so-called 
UV-renormalons In QCD, it is the low momentum integration region of virtual 

gluons dressed with running couplings which leads to non integrable singularities in the Borel 
plane the so-called IR-renormalons |3|,|6|,|l3fl. 

The study of renormalon properties in gauge theories is at present an active field of 
research. The role of IR-renormalons in the operator product expansion of two point functions 
and their relationship with non-perturbative inverse power corrections has been extensively 
discussed in the literature f9Hi~6| . This and further discussions which originated after ref. |17| 
appeared, has led to a new point of view [^] concerning renormalons in QCD which focuses on 
the possibility that their systematic study in a given hadronic process might suggest generic 
non-perturbative effects of a universal nature. This applies to the case of IR-renormalons as 
well as to the much less explored role of UV-renormalons [19, 2(J in QCD. 

Practically all work on renormalons in the literature is restricted to the effect of the 
first coefficient f3\ of the /^-function. There is a good reason for that : the position of the 
singularities for positive Borel variable is governed by (3\\ and it is the sign of (3\ which fully 
dictates their fate as UV-like {(3\ > 0) versus IR-like singularities {(3\ < 0) . It is also often 
stated that these singularities are due to the Landau pole in the running coupling constant. 
The purpose of this note is to discuss some interesting properties which appear when the first 
two terms of the /3-function are fully taken into account, and to show how they clash with 
this common belief that renormalon singularities are due to the Landau pole in the running 
coupling constant. 



2. We shall first review the usual analysis of renormalon effects in the precise case of 
the anomalous magnetic moment of the electron a e in QED, the same example that Lautrup 
considered in ref. The relevant Feynman diagrams are those generated by the one- 
renormalon chain in Fig. 1, and their contribution to a e is given by the integral 



(1) 



where x is the Feynman parameter which combines the propagators in the loop vertex of 
Fig. 1, and a e Q the QED effective charge. For large values of the euclidean momentum 
k 2 (k 2 /m 2 = in our case, with m the electron mass) the effective charge obeys the 
renormalization group equation : 



^For a comprehensive review of the subject and a collection of early articles see ref. |^| 
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where the infinity superscript refers to this asymptotic limit for k 2 . 

Equation is an exact integral representation of the contribution to the electron anomaly 
from the infinite class of diagrams which Fig. 1 represents. In particular, if one takes a e g = 
a ~ 1/137 one of course re-obtains the value a e = a/2ir. When the number of vacuum 
polarization bubbles is large, the integral is then dominated by the x — * 1 region; i.e. large 
values of the virtual euclidean loop momentum, and becomes equivalent to the integral 



1 r°° dk 2 ( m 2 \ 2 



vr Jm 2 k 2 \ k 2 




(3) 



The insertion in the integrand above of the one loop solution to the renormalization group 
equation in (|2|) with (3{a) = = |]f ; i.e., the insertion of the expression 

a Pi k 2 

l a log (4) 



'eff (£' a ) 
in (||) , and the change of variables : 



a (°°) ( ML n \ 2ir 6 m 2 



B- = 1 -, — , , r i where B = — , (5) 
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leads then to the Borel sum 



oo 



dze a — — . (6) 

vr Jo 2-Bz v; 

The singularity (UV renormalon) at z = 2/B reminds us that we only know the function 

for z < 2/B, just as we only know (j^z,a(rn)^ in Eq. (Q) up to a momentum k 2 = 

A| = exp (j^j ■ At this momentum (the Landau scale) the effective charge becomes infinity. 
The association of the UV-renormalon ambiguity with the Landau pole comes from this 
observation. 

Notice that in Eq. ^ the leading UV renormalon at z = 1/B does not contribute, so 
that the ambiguity starts at z > 2/B. The reason for it is the factor 1 — x in the numerator 
in eq. (|l]), which appears because of the electron helicity conservation of the electromagnetic 
interaction, and leads to 2-powers of the factor m 2 /k 2 in the euclidean momentum version 
of the same integral in (||). In the language of UV-effective operators of Parisi 0, this 
corresponds to the tree level insertion of the operator (m 3 /A| / ) tj;(x)a f j iU ip(x)F^ L ' (x). 



3. Let us now study the same a e observable but in the presence of the first two terms of 
the QED /^-function : 

/3H = [/3i = |]- + [/?2 = M-) 2 . (7) 

In this equation we have no loss of generality since any (3— function can be brought into this 
form by a (perturbative) coupling constant redefinition. 
To simplify the notation we shall denote 
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The general solution of eq. (pi) is then 



a(m) . . k 
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Ba(m) log — g + 8Ba(m 



a{k) m 



a(m) 1 + 5Ba(k) 



log — 77T + 1°§ 



a(k) 1 + 5Ba(m) 



(9) 



where a(m) is the boundary condition at k = m . The shape of these solutions in the plane 

l°g"^r> oi(k) is shown in Fig. 2. In full generality, and depending on whether the boundary 

condition is chosen to be a{m) > 0; —jg < a(m) < 0; or a(m) < — the solutions lie 

along the curves shown in the regions I, II, or III of Fig. 2. These regions correspond to the 

analogous ones in the /3-function of Fig. 3 where the arrows show the flow of the effective 

charge as k 2 — > oo. The three regions are separated by the fixed points at a* = and 
** _ _j_ 

a — SB' 

The physical region is clearly the one in I. For a given initial condition a{m) > 0, the 
running coupling a(k) grows and hits a Landau pole at a finite value of k 2 that we call A^, 
where a(A^) — > oo . With an initial condition fixed in region I, there is no solution to the 
renormalization group equation when k 2 > A 2 . Yet the integral over the euclidean momentum 
in eq. (0) generated by perturbation theory runs up to oo. It has been recently shown by 



Grunberg [21| that in the presence of the first two terms in the power expansion of the f3— 
function it is still possible to find an exact change of variables which maps k 2 integrals in the 
euclidean space, like the one in (||), to integrals in the z-variable of the Borel plane. The 
change of variables in question is 



1 



a(m) 



B \ = 1 + SB^m) ' (10) 



and leads to the Borel sum 



1 f°° e~ zSB 

a e = — / dze " s — - — r. (11) 

2W (1-^Y+ S V ' 

However this is only a formal expression. With the initial condition a{m) > fixed, the 
change of variables in ( |l0|) is only well defined in the region where m 2 < k 2 < A| , which 
corresponds to < -B| < (1 + <5i?a(m)) _1 < 1 . Let us however inspect what happens 
when one goes beyond this physically allowed region i.e., when (1 + SBaim))^ 1 < i?| < oo. 
This is equivalent to taking a certain arbitrary a(k) in the non-physical region III, where 
(1 + 5Ba(m))~ l < B^ < 1, and then in the region II, where 1 < -B| < oo . With this proviso, 



cq. (11) has a definite (but arbitrary) meaning. From this point of view, the singularity at 
z = 2/B acts as a reminder of this arbitrariness. For example, one can define an "extended" 
a(k) valid at all fc 2 -values, m 2 < k 2 < oo, such that in the physical region m 2 < k 2 < A| it 
coincides with the physical a(k) and that it satisfies the perturbative renormalization group 
equation in (|9|), except at a finite number of singularities. This extended a(k) has the following 
forrr|]: 



aim) _ „ , . , k 
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1 — Ba(m) log — j + 8Ba(m) log 



a(m) (l + 5Ba(k)) 



a(k) (1 + 5Ba(m)) 



(12) 



a(k) v ' ° m 2 

This expression leads to the trajectory with the arrows displayed on Fig. 2, where the arrows 
correspond to the flow as one integrates over z going first from region I to region III, and 



'"Notice the absolute value in the argument of the logarithm, in contrast to eq. 
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then to region II. In so doing, one sees that one crosses three singularities : one in which 
a(k) goes to infinity -the conventionally called Landau singularity- and two in which log 
goes to infinity -at the fixed points of f3(a) : a* = and a** = — . With this "extended" 
definition of ct(k), the change of variables in ([n]) leads to the Borel sum 



I roo „—z6B 



dze q j. (13) 

2nJ 



:i-f: 



The two loop Borel integrals (11) and ( |l3| ) still have the singularity at z = 2/B, as in the one- 
loop case. This may lead one to believe that, exactly as in the one-loop case, this singularity 
is caused by the corresponding Landau singularity in euclidean momentum. However this is 



not so. To see this explicitly let us expand eq. (||), or eq. fll2j), for a(k) — > +00 i.e. near the 
Landau pole : 

k 2 { 1 \ a 1 / 1 \ 

B a lo S - 2 = 1 - Htalog (l + m ) - + O (-^) . (14) 

The terms of O f ^gy) cancel, which means that the Landau pole of a(k) at 

t . 2 = Ai = m2exp {_L_ 5l0g ( 1 + _i_)}, (15) 

goes like 

a(k) = ^ (io g A|-iogP)V2 ; ( 16 ) 

i.e., it is a squared root singularity and therefore integrable up to the Landau scale where 
Bz/2 = (l + 5Ba(m))~ l < 1. From that point onwards the integration in z can be done using 
the extended a(k) in eq. (12), following the arrows on Fig. 2 from region I to region III until 
a = a** where one hits the Borel singularity at z = z n . This is to be contrasted with the one 
loop behaviour in eq. (|j) where 

° (A;) = llogA|-log^' 

and the singularity is not integrable. The euclidean origin of the z = 2/B singularity in the 
Borel sum at the two-loop level is the log singularity corresponding to a** = — ; not 
the Landau pole ! Of course in the limit /?2 — ^ 1^ one recovers the one-loop situation in which 
a** —> 00, and the Landau pole. 



4. The basic features we have discussed in the case of the anomalous magnetic moment of 
the electron in QED are rather generic and they appear as well in the case of IR-renormalon 
calculus in QCD where, typically, one encounters euclidean integrals like 

rQ 2 rlk 2 f k 2 \ n / k \ 

*-w>=j - (!■•«•«»)■ < is > 

§ Although one can invoke 1/Nf arguments to effect the limit /J2 — » in QED, there is no analogue in the 
case of QCD. 
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with n > 0, so that the integral is infrared convergent. The scale Q 2 in these integrals 
corresponds to a sufficiently large choice of euclidean momentum at which the QCD running 
coupling a s (Q) is reasonably small. Here the ambiguity problem appears because of the 
integration over virtual euclidean momentum in the infrared region < k 2 < Q 2 , where the 
extrapolation of the perturbative a s (Q) coupling is not well defined. It is often stated in the 
literature that the reason for the ambiguity is the existence of the Landau pole in the region 
of integration. We shall see that, as already shown in the case of a e in QED, this is only 
correct if one rather arbitrarily restricts the QCD beta function to its first term. 

The discussion runs parallel to the one in the previous section, except that now the running 
coupling a s a s (Q)j = a s (k) obeys the QCD renormalization group equation 

^^ = -b a 2 M-bial(k), (19) 
dlog|p 

where we are using the notation 



^-^/3! QCD = ^(HN c -2n f ), (20) 

i _ 1 n QCD 1 J 17 2 N 2 -l 5 \ 

01 = k Po = — k \ — N„ nt ±y c nf ) ■ 21 

2vr 2 H2 8vr 2 \ 3 c 2N C 1 3 ; J v 1 

with Do, &i > in the physically interesting case N c = rif = 3. We can consider eq. (|l9|) to 
be a generic case since any other /3-function can be brought into this form by an appropriate 
coupling constant redefinition. Of course this redefinition will entail in general an infinite 
power series. The solutions of eq. ([H]) have the pattern shown in Fig. 4, with region I 
corresponding to the physical region with the boundary condition a s (Q) > 0. The equivalent 
change of variables to the one in eq. ( Jl0| ) is [21| 



I _ <Xs{Q) 

Zn~ l + |a s (Q) ! Zn = T - 
It maps the k 2 integral R n onto an integral in the z-variable of the Borel plane 



(22) 



MQ 2 ) = --l «. = ^. (23) 



For a fixed n, the Borel integral ( |23[) is singular at z = z n , which according to ([22|) happens 
when a s (k) = — ^ and has nothing to do with the Landau pole at which, by definition, 
a s (k) = oo . In fact at the Landau pole, which now occurs when 

k 2 = A£ = Q 2 exp(-4- + felog(l + — ^-H , (24) 



\b a s (Q) ^ b 2 S V ^ oia s (Q) 
a s (k) — » oo as a squared root singularity : 



° s(A:) " (260V2 (log A| _ bg jfe2)i/2 ' ( 25 ) 



and it is therefore integrable. 
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Notice that the ambiguity in eq. (p3| ) can be parameterized by a power-like contribution 
in Q 2 . Indeed, making the change of variables 



± = 1 + aiQ h )UJ (26) 

one finds that the ambiguity in R n (Q 2 ), 5R n (Q 2 ), is given by 

SMQ) J^i e -^ f^U ' 4)\-* (27) 



n J* uj 



J \a(Q) b 



/AM 

5R n {Q 2 ) = {coefficient } x ( -± J (28) 

where / t cL; e _a; denotes an integral with an arbitrary prescription to skip the sin- 

gularity at u = 0. Although, as repeatedly emphasized, the Landau singularity is not directly 
responsible for the singularity at z = z n , the final result of eq. (^) is insensitive to this fact 
and the form of ( |2"8| ) is the same as in the case of the one-loop (3 function. 

Before we conclude, we would like just to mention that an analogous reasoning for a e in 
QED leads, mutatis mutandis, to an equation like ( |28| ) for the analogous quantity 5a e but with 
n = —2 and Q 2 = m 2 . Therefore it is also true for a two-loop (3 function that the ambiguity 
Sa e can be given by the tree-level insertion of the local operator (m 3 /A| y ) ^(x)a fJiU ip(x)F fl ' / (x), 
just as in the one-loop case. 
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Note Added: 

The analysis of renormalon effects in the case of a j3 function with a perturbative fixed 
point has been considered by Grunberg in ref. [21] and also by Yu.L. Dokshitzer and N.G. 
Uraltsev in hep-ph/9512407 . They reach similar conslusions to ours. We thank N.G. Uraltsev 
for drawing his paper to our attention. 
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Figure Captions 



Fig. 1 : Feynman diagrams which define the one-renormalon chain contribution to a e in 
eq. ©. 

Fig. 2: General pattern of the QED renormalization group solution in eq. @. r stands for 
log-^r. The regions I, II and III correspond to the choice a{m) > 0; — < a(m) < 0; and 
a(m) < — as boundary condition. 

Fig. 3 : The QED /^-function in perturbation theory at the two-loop level. The arrows show 
the flow of the effective charge as k 2 — > oo. 

Fig. 4 : General pattern of the QCD renormalization group solution in eq. (|l9|). r stands for 
log^-. The regions I, II and III correspond to the choice a s (Q) > 0; — j& < a s (Q) < 0; and 
&s(Q) < — ^ as boundary condition. 
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